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ABSTRACT 



In the present work we show robust indications of the ex- 
istence of g modes in the Sun using 10 years of GOLF 
data. The present analysis is based on the exploitation of 
the collective properties of the predicted low-frequency 
(25 to 140 /iHz) g modes: their asymptotic nature, which 
implies a quasi equidistant separation of their periods for 
a given angular degree (£). The Power Spectrum (PS) of 
the Power Spectrum Density (PSD), reveals a significant 
structure indicating the presence of features (peaks) in the 
PSD with near equidistant periods corresponding to £=l 
modes in the range n=-4 to n=-26. The study of its sta- 
tistical significance of this feature was fully undertaken 
and complemented with Monte Carlo simulations. This 
structure has a confidence level better than 99.86 % not 
to be due to pure noise. Furthermore, a detailed study 
of this structure suggests that the gravity modes have a 
much more complex structure than the one initially ex- 
pected (line-widths, magnetic splittings...). Compared to 
the latest solar models, the obtained results tend to fa- 
vor a solar core rotating significantly faster than the rest 
of the radiative zone. In the framework of the Phoebus 
group, we have also applied the same methodology to 
other helioseismology instruments on board SoHO and 
ground based networks. 

Key words: Sun: Helioseismology, Sun: Interior, Sun: 
Oscillations. 



tic (p) eigenmodes and the subsequent structural inver- 
sions, providing the stratification of crucial variables like 
the sound speed down to 0.05 Rq. However, the acous- 
tic modes are less sensitive to other structural variables 
such as the density. In this case, the agreement with the 
models is still poor in the deepest layers of the radiative 
zone iH. Moreover the dynamical properties of the so- 
lar interior are practically unknown in the nuclear core. 
Indeed, due to the limitation in the number of acoustic 
modes penetrating this region and to their low sensitivity 
there, the solar rotation profile (see Fig.[T]i is very uncer- 
tain below 0.2 Rq [e.g. 2, 3, 4]. At these depths, large 
uncertainties also exist because of their lack of sensitiv- 
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1. INTRODUCTION 

Detailed knowledge of the solar interior has been ob- 
tained by means of precise measurements of its acous- 



Figure 1. 2-D inversion of the solar rotation profile using 
the modes t < 25 from GOLF [3] and MDI [5] instru- 
ments. The horizontal error bars increase towards the 
centre of the Sun. Below 0.2 the rotation profile is un- 
known. 



ity. On the contrary, solar gravity (g) modes exclusively 
provide information on the radiative region and the inner 
core. They are required to obtain a complete picture of 
the interior of our Sun. Unfortunately, solar g modes have 
not yet been clearly detected although extensive and in- 
tensive searches have been developed along past decades. 
These modes are evanescent in the convective zone and 
their amplitudes at the solar surface could be very small 
To look for the g modes, two complementary ways 
of searching strategies have been developed during the 
last decades: 



• The first method consists in looking for individual 
modes (relevant spikes in the Fourier spectrum of 
the observed signal) above a given statistical thresh- 
old (typically 90 % confidence level) at frequencies 
above 150 /iHz. The first studies from the GOLF 
(Global Oscillations at Low Frequency) Team |@] 
and the Phoebus |9] and BiSON (Birmingham Solar 
Oscillations Network, [10]) groups looked for indi- 
vidual peaks (spikes). An upper limit of 1 cm/s has 
been found for these single peaks. Another study of 
the GOLF data was devoted to the search for multi- 
plets - instead of individual peaks - which has low- 
ered the detection level to more than a half. Two 
candidates have been observed and followed with 
time with more than 90% confidence level in a 20 
yuHz interval [11] and one of them has even been de- 
tected with a 98% confidence level after 3000 days 
of observations in a 10 /iHz interval fl^. However, 
some ambiguity remains on the proper labelling of 
the detected components, in term of their £, m and n 
values. Nevertheless, some scenarios have been pro- 
posed in order to explain the visible peaks which can 
constrain the physics and dynamics of the solar core. 
Indeed, the frequencies of such modes are extremely 
sensitive to the structural properties of the radiative 
zone, including the problem of the solar abundances 



• An alternative way, developed in the early 80s, is 
to look for t he g -mode asymptotic properties (see 
for example ill4i ITsll and references there in): this 
theory predicts a constant separation (in period) be- 
tween the central components of the gravity modes 
at low frequency. Fig.|2]shows the separation in pe- 
riod of consecutive g modes (APf = Pi^{n+i) — 
Pi,n), with £ < 3 and periods below 11 hours, 
as a function of the period computed from a seis- 
mic solar model [2] assuming a rigid core rotation 
(Wa . The effect of the rotation is to split the modes 
into different m-components regularly spaced in fre- 
quency. As a consequence, only zonal ones, m=0, 
are equidistant in period while the others shift pro- 
gressively as the period of the modes increases. 
Consequently, this method becomes extremely sen- 
sitive to the rotation rate of the inner solar layers. 
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Figure 2. Separation in period of theoretical g modes 
with a rigid rotation rate in the core for £=1 (red), £=2 
(green) and 1=3 (blue) modes. The red stars correspond 
also to the same £=1 separations but for a model with a 
rotation rate in the core 5 times larger than in the rest 
of the radiative region. The region inside the vertical 
dashed lines corresponds to the region of interest (25 to 
140 fiHz). 



2. ANALYSIS TECHNIQUE 



Using the g modes showed in Fig. |2] and considering that 
they have long enough (more than ten years) lifetime and 
a constant amplitude for each mode degree (with the same 
relative visibility ratio than for the p modes), we can build 
up a simple model - without noise in a first step- of the 
power spectrum between 25 and 140 /iHz, PSD{v) (see 
top panels of Fig. |3]l. To analyze the existing periodici- 
ties inside this spectral region, the independent variable 
of PSD{v) is changed from frequency to period, giv- 
ing PSD{P) (middle panels of Fig. |3]l. Then, the pe- 
riodogram of the PSD{P), hereafter called PS, is com- 
puted using a 10-time-oversampled sine- wave fitting al- 
gorithm (bottom panels of Fig.[3]i. The set of modes cor- 
responding to each degree gives rise to a peak on the PS 
and to a sequence of harmonics of it. For the £=1 modes 
the main peak in the PS is centered at APi ~ 23.5 min- 
utes with an asymmetric shape because the separation of 
the modes at lower periods decreases (see also Fig. |2|. 
The behavior of the £=2 and 3 modes can be explained in 
the same way. Looking at the PS including all the modes 
together (right bottom panel of Fig. [3]), only the APi can 
be easily discerned. The others, AP2,3, are mixed with 
the harmonics of the £=l and it is not easy to recover 
them. The presence of the £=2 modes have also an influ- 
ence on the APi due to the beating between them. The 
contribution of the £ — 3 modes is negligible due to their 
small simulated visibility. For these reasons we have de- 
cided to concentrate our study only on the signature of 
the APi between 22 and 26 minutes. It is also important 
to notice that a change in the relative amplitude between 
the different £ modes will modify the relative amplitudes 
of the peaks of each AP^ . 
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Figure 3. Simulated PSD of theoretical g modes computed from the seismic model using a rigid core rotation law for i=l 
modes (left panels) and I < i < 3 (right panels) expressed in frequency (top) and in period (middle). The bottom panels 
correspond to the PS of each above spectrum. In the bottom left panel a clear peak appears centered at 23.5 minutes 
corresponding to the main periodicity of the £ ^ 1 modes as well as all the harmonics at lower periods. In the bottom 
right panel the structure of the main separation of the 1=2 modes appears centered at 14 minutes. 
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Figure 4. Top: Zoom between 2 and 8.5 hours of the 
reconstructed waves in the PSD(P) using the g modes 
computed from the seismic model with a rigid core ro- 
tation. BottomiSame as before but using a core rotation 
rate 5 time faster than the rest of the radiative region at 
R = 0.15Rn. 



From the APi pattern measured in the PS, it is possible to 
reconstruct the waves in the PSD(P) corresponding to the 
peaks between 22 and 26 minutes and its first harmonic. 



Thus, we compute: 

A,sin(^ + 0J+ ^ A,sm(^ + (j) 

Pie[22:26] ^ * ^ P,e[iia3] ^ ' 

(1) 

By doing so, we are also taking into account the phase 
of the waves in the PS. This information constrains the 
position of the maximum of the sinusoidal waves that 
should coincide with the positions of the g modes in the 
PSD. Fig.Hlshows the reconstructed waves superimposed 
on the original £=l simulated modes for the same solar 
seismic model but with two different core rotation laws. 
Firstly, the sinusoidal waves reproduce the pattern of the 
g modes and they are very sensitive to the splitting of 
the modes (see the bottom panel) where the waves are 
divided in two at 5.5 hours to follow both split compo- 
nents. Thus, this method can be directly used to check the 
general trend of the solar rotation profile inside the core. 



3. DATA ANALYSIS AND DISCUSSION 



Data from the GOLF instrument aboard the SoHO mis- 
sion are used in this work. GOLF is a resonant scattering 
spectrophotometer measuring the line-of-sight Doppler 
velocity of the sodium doublet |17]. We use 10 years 
time-series starting 1996 April 1 1 and properly calibrated 
into velocity following II18I1 . After computing the corre- 
sponding Fourier Transform, the solar background was 
properly modeled and removed by means of a first order 
polynomial fit in a log-log scale. The resultant PSD{v) 
is flat in the region of interest and its statistical properties 
are well represented by a with 2 degrees of freedom 
(d.o.f.). The resultant PS is plotted in the top panel of 
Fig. |5] Between 22 and 26 minutes a wide pattern of 
high amplitude peaks appears. The maximum amplitude 



in this region reaches 6.5-ct while the average power is 
2.95 times higher than the averaged power of the rest of 
the spectrum. Comparing this pattern with the one pro- 
duced by the theoretical £=l modes issued from the seis- 
mic model (bottom panel of Fig. |5]l we find that the PS 
are both alike. 
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Figure 5. Top: PS of the PSD expressed in period com- 
puted from the GOLF velocity time series. Bottom: PS 
of the PSD{P) from theoretical g-modes computed from 
the seismic model and using a rigid core rotation. 

We have studied the PSD{v) between 25 and 140 /xHz, 
because this region is a compromise between the need 
of having modes in the asymptotic region and with high 
enough amplitudes (at lower periods, see: f?', T^). We 
have verified that when shorter spectral regions are used 
in the calculation of the PS, e.g., from 30 or 35 to 140 
/iHz, the prominent structure around 22 to 26 minutes 
still appears in the PS with maxima at 7.3- and 6-cr re- 
spectively (see the three panels of Fig.|6l). 

To check the probability that this kind of structure placed 
between 22 and 26 minutes (guided research) is due to 
pure noise (Ho hypothesis test) we compute a Monte 
Carlo simulation of 6 10^ realizations (uncertainty on the 
results of 0. 13 %) of a pure noise power spectrum follow- 
ing a with 2 d.o.f. ifTotl . To characterize this feature, 
we use two indicators: the maximum power in the region 
between 22-26 minutes (6.5-<j) and the amount of total 
power in this region (2.95 above the rest of the spectrum). 
Thus, with 99.49 % of confidence level, we reject the i/o 
hypothesis. In the framework of the PHOEBUS group 
II20I1 we also compute the probability to find this kind of 
pattern but anywhere in the PS (non guided research) be- 
tween 5 and 50 minutes. In this case, with a confidence 
level of 94.15%, the Ho hypothesis is rejected. 

The Monte Carlo simulations show that pure noise is very 
unlikely to generate this kind of patterns, then we can 
try to check if they can be due to periodic signals in the 
PSD placed where we expect the gravity modes to be. 
Therefore, we reconstruct the waves in the PSD(P) corre- 
sponding to the APi and its first harmonic. The result is 
plotted in Fig.ll] To characterize it we compute the cor- 
relation of the resultant waves with the ones computed 
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Figure 6. PS of three different regions of the PSD, from 
top to bottom 25, 30 and 35 fj,Hz to 140 iiHz respec- 
tively. The cyan shaded zones correspond to where the 
APi structure is expected, between 22 to 26 min, and its 
first two harmonics. The orange shaded zones are the 
same but for the m = component of the AP2. 
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Figure 7. Reconstructed waves using GOLF data com- 
pared to the reconstructed waves from the g modes of the 
seismic model using a rigid core rotation law ( top) or us- 
ing a 5 times faster core rotation law (bottom). 



from solar models changing several parameters: the solar 
model, the rotation law and the solar core rotation axis 
inclination. In most of the cases we see that, for a given 
solar model, the correlation increases with a higher rota- 
tion rate. It reaches a maximum between 3 to 5 times the 
rotation of the rest of the radiative region, depending on 



the other parameters taken into account. This can be seen 
by comparing both panels of Fig. |7] In the top one, the 
GOLF reconstructed waves start to be larger at 3.8 hours 
and are clearly split in two different waves at 5.5 hours 
while the model with a rigid rotation has single thinner 
peaks all over the plotted region. In the bottom panel the 
correlation between both curves is higher and their gen- 
eral behaviours are similar. Only above 7.5 hours the 
GOLF waves shift to higher periods in comparison to the 
simulated ones. 

Using the same Monte Carlo simulation as before, we can 
check again the Hq hypothesis comparing the resultant 
reconstructed waves with a model an imposing a thresh- 
old of only 20% which is the smallest one we found when 
comparing the real GOLF data with the models. In this 
case Hq is rejected with more than 99.86% confidence 
level. 

Up to now, we have seen that the studied structure is un- 
likely due to pure noise. One could think that it may be 
due to some kind of instrumental noise. However there 
are no known instrumental effects having periodicities in 
the range of 2-14 hours and being equidistant in period. 
The other possibility is that these periodicities have a so- 
lar origin. Excepting a possible relation with convection 
we can analyze the consequences of assuming that they 
are produced by g-mode asymptotic properties. With 
such objective, we have introduced noise - a with 2 
d.o.f. - to the simulations used previously. The resultant 
PSD (in arbitrary units) and corresponding PS are shown 
in Fig.lSt and Fig.|8}; respectively (for a simulation with 
a Signal-to-Noise ratio (S/N) of 30. The visibility of the 
APi structure is degraded because of the noise. Instead 
of having a big peak - like in the noise-free case - we 
have a more complicated structure (the exact distribution 
of peaks depends on each noise realization). However, 
in order to keep a significant APi structure (above the 
limits defined in the Monte Carlo simulations) we need a 
very high signal/noise (S/N) ratio in the power spectrum, 
much higher than the GOLF observed one which seems 
to be around one. Otherwise we should have already mea- 
sured the g modes. 

Therefore, it is not possible to find, in average, a proper 
APi - with a small S/N ratio in the power spectra - 
unless the probability to have energy on the modes is 
transfered into more bins. We have found three ways to 
reach such objective: a) assuming a finite lifetime of the 
modes as suggested by the recent works by |,21,] which 
spreads the power of each mode from a single bin into 
a Lorentzian profile covering a larger region; b) increas- 
ing the visible multiplet components assuming a different 
rotation axis inclination of the core than in the radiative 
region, as it was first suggested by c) considering 
a magnetic field in the core or at the base of the convec- 
tive zone. Indeed, |23] and more recent calculations fv^ 
have given a complicated picture of the contribution of 
a deep magnetic field with many components appearing 
for each mode. To explore these possibilities we have 
chosen the first. Thus, Fig. [8J5 corresponds to the power 
spectrum of the same simulation as in Fig. [8^ but using a 
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Figure 8. PSD of two different simulations (a, b) and 
their corresponding PS spectra ( c, d). See the text for 
more details. 

finite lifetime - several months - of the g modes. In this 
case the S/N can be significantly reduced (S/N=5) while 
the PS continue to show a detectable APi structure (See 
Fig. |8}l). To go further, a more statistical approach of 
these simulations is needed including heavy Monte Carlo 
simulations covering the different physical mechanisms 
explained above. 



4. APPLICATION TO THE VIRGO/SPM DATA 

In the framework of the PHOEBUS collaboration H^], 
we have applied the same methodology to 3473 days of 
the sum of the 3 SPM channels of the VIRGO experi- 
ment. Unfortunately, no clear structures appear in the PS. 
This can be a consequence of the different solar back- 
ground comparing disk-integrated velocity and intensity 
data. 

In order to go a step further, we put to zero all the points 
above or below a given threshold before computing the 
PS. Left panel of Fig. |9] shows the maximum of the high- 
est peak between 22 and 26 in the PS when the level of the 
threshold is continuously increased (or decreased). For 
example, the maximum is between 2 and ^ 4 when we 
put to zero all the points above a threshold changing from 
2 to 9 (continuous line in the figure). However, when we 
remove those peaks in the PS with lower power (less than 
the threshold, dashed line in the figure) we see that for 
a threshold between 5.5 and 6.5, we have a peak in the 
studied region with an amplitude of ~ 6. Right panel 
of Fig.|9]shows the PS computed from a PSD where only 
peaks with S/N ratios within the range 5.5 to 10 have been 
considered. A similar peak-structure to those found us- 
ing GOLF data can be clearly seen between 22-26 min. 
The pattern seems to be larger than the GOLF one. As 
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Figure 9. Left: Power of the highest peak between 22 and 26 minutes in the PS when the VIRGO/SPM PSD is thresholded 
above (continuous line) or below (dashed line) a given value between 2 and 9. Right: PS corresponding of the PSD 
computed using 3473 VIRGO/SPM intensity data. A minimum and a maximum threshold have been applied to the peaks 
in the PSD between 5.5 and lO-a. See the text for details. 



a consequence of the threshold done, only 50 bins in the 
PSD are responsible for this structure (the PSD of the 
Virgo data between 25 and 140 /iHz has more than 32000 
points). Thus, a visual identification of these peaks can 
be directly done in the PSD. 

Further studies should be carried out by the PHOEBUS 
group before fully understand the VIRGO/SPM data. If 
the signal is confirmed, it rules out a GOLF instrumental 
origin. However, as both instruments are on the same 
spatial platform an external confirmation (i.e. BiSON or 
GONG data) would be extremely useful. Finally, the deep 
study of the non disk-integrated MDI/SoHO data would 
bring the opportunity to access higher degree modes with 
new valuable information on the inner core. 



5. CONCLUSIONS 

In the present work we have shown the presence of pe- 
riodic signals in the PSD of the GOLF spectrum when 
it is expressed in period. We have shown by the Monte 
Carlo simulations that it is very unlikely that these pat- 
terns were due to pure noise with 99.49 % confidence 
level. Assuming that it is due to g-mode signals, it can 
only be produced by the presence of £=\ modes with 
—4 < n < —26. A direct comparison with theoretical 
g modes showed that the Hq hypothesis can be rejected 
with more than 99.86% confidence level. 

From this analysis we start to have access to the dynam- 
ics of the internal core and we have shown that the core 
seems to rotate in average between 3 to 5 times faster than 
the rest of the radiative zone when unrealistic "rotation- 
step" profile is used. More work should be done in the 
future to introduce more realistic core rotation law pro- 
files and also to use the information corresponding to the 
higher degree modes. 



Another important consequence of the analysis per- 
formed here is that the individual structure of the g-mode 
multiplets should be more complicated than what, up to 
now, has been thought. To simulate realistic PSD (low 
S/N) we need a higher number of bins containing signal. 
As we have seen, at least, three different scenarios are 
compatible with a significant reduction of the required 
S/N. Further studies will be necessary to discriminate be- 
tween them and to improve our knowledge of this modes. 
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